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Abstract. In atomic force microscopy (AFM) tip-surface interactions are usually 
considered as functions of the tip position only, so-called force curves. However, 
tip-surface interactions often depend on the tip velocity and the past tip trajectory. 
Here, we introduce a compact and general description of these interactions 
appropriate to dynamic AFM where the measurement of force is restricted to 
a narrow frequency band. We represent the tip-surface interaction in terms 
of a force disk in the phase space of position and velocity. Determination of 
the amplitude dependence of tip-surface forces at a fixed static probe height 
allows for a comprehensive treatment of conservative and dissipative interactions. 
We illuminate the fundamental limitations of force reconstruction with narrow 
band dynamic AFM and we show how the amplitude dependence of the Fourier 
component of the force at the tip oscillation frequency, gives qualitative insight into 
the detailed nature of the tip-surface interaction. With minimal assumptions this 
amplitude dependence force spectroscopy allows for a quantitative reconstruction of 
the effective conservative tip-surface force as well as a position-dependent damping 
factor. We demonstrate this reconstruction on simulated intermodulation AFM 
data. 
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1. Introduction 

For the understanding of surface reactions and the characterization of materials it is 
desirable to measure local forces close to a sample surface. The most common method 
to measure these surface forces is atomic force microscopy (AFM) [lj . Historically, the 
first force measurements were static measurements for which the force is presented 
as a scalar function of the static tip-sample separation, the so-called force curve[2j E]- 
This representation is sufficient for conservative forces but the total tip-surface 
force may also contain contributions from dissipative forces. Since dissipative forces 
depend on probe velocity and past trajectory, dynamic force spectroscopy methods 
are required for their measurement. Moreover, the visualization of dissipative forces 
as a function of position is valid only for a specific probe trajectory and simple force 
curves cannot capture the full character of the interaction. Despite the development 
of several dynamic methods jU EJ El EJ El [TOj [H] surface forces are still usually 
treated as functions of the probe position only and represented by simple force curves. 

Here, we present a comprehensive framework for the representation and analysis 
of complex surface forces as they are measured by dynamic AFM. We concentrate 
on the most common modes of dynamic AFM: amplitude-modulated AFM (AM- 
AFM) and frequency-modulated AFM (FM-AFM), which can be considered as 
narrow frequency band methods|12j. We explore the fundamental limit of force 
reconstruction with narrow band dynamic AFM at fixed probe height and show 
how minimal assumptions allow for a quantitative reconstruction of the tip-surface 
interaction. 

1.1. Cantilever-based dynamic force measurements 

At the heart of the AFM apparatus is a micro-cantilever with a sharp tip. The 
cantilever is firmly clamped at one end and the tip is located at the other end which 
can move freely. It is assumed that surface forces only act on the tip whereas the rest 
of the cantilever does not experience significant surface forces. In dynamic AFM an 
additional external drive force is applied to maintain an oscillatory motion. Thus, 
the dynamics are governed by the force between tip and surface, the external drive 
force and the properties of the cantilever beam. 

Since the cantilever is a three dimensional continuum object its motion is usually 
described by the amplitudes of different oscillation eigenmodes. In general, these 
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modes can cause the cantilever to bend in all directions in space. However, the 
cantilever is positioned such that the softest flexural modes bend the beam in a plane 
orthogonal to the surface plane. We restrict ourselves to the case where only these 
flexural modes are excited by the drive force. Due to this experimental configuration 
the cantilever is much more susceptible to the component of the tip-surface force 
which is orthogonal to the surface plane. This component of the force is typically the 
most dominant component and the influence of lateral force components is considered 
negligible. In this case the cantilever acts as a mechanical projector which reacts only 
to one component of a three dimensional force vector field. 

The deflection w of a cantilever of length L orthogonal to surface is described 
by a one dimensional Euler-Bernoulli equation |13] 

EI ?^ + ^^A = F{xM (1) 

where E is the Young's modulus, I is the second moment of area, fi is the mass per 
unit length of the cantilever, x is the position coordinate along the cantilever beam 
and t is the time variable. The force term F includes the surface forces acting as 
a point-like load at position x = L, the external drive force and the hydrodynamic 
damping due to the surrounding medium [13] . To express the solution to equation 
([1]) in terms of the oscillation eigenmodes, one performs a separation of variables and 
obtains a solution for w(x,t) that is a linear combination of different mode shapes 
$( n )(x) with time-dependent amplitudes q^ n \t), 

oo 

w (x,t) = J2q {n) (t)® {n \x) (2) 

n=l 

where the mathematically orthogonal modes $ ( - n ^ have different resonance frequencies 

(n) 

In most dynamic AFM modes only the first eigenmode is externally excited. 
The spectrum of the resulting tip motion is then confined in a narrow frequency 
band around the first flexural resonance frequency Wo^^J. Due to the fact that 
the resonance frequencies of higher eigenmodes are not integer multiples of the first 
resonance frequency and that the eigenmodes have high quality factors, only the 
first mode contributes to the cantilever motion and the tip motion z(t) = w(x = 
L,t) + h can be approximated as 

z(t) + (3) 
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where h is the static probe height above the surface. The time-dependence is given 
by an effective harmonic oscillator equation [151 CE] 



where the dot denotes differentiation with respect to the time t, is the quality 
factor of the first flexural resonance, the effective mode stiffness, -Fdrive the time- 
dependent external drive force and F ts the force between tip and surface. Thus, 
equation (00) reduces the basic physics of dynamic AFM to that of a damped harmonic 
oscillator which is moving in an nonlinear force field F ts and is subject to a time- 
dependent drive force i*d r i ve . 

Simulations starting from a given tip-surface force reveal multiple oscillation 
statesp/fl [18] and period multiplications [19], qualitative features of the dynamics 
which have also been observed in experiments |20l |21~1 122]. However, the fundamental 
challenge in AFM is actually the inverse problem: Given an accurate and well- 
calibrated measurement of the dynamics z(t), how can we quantitatively determine 
the tip-surface force. Different sophisticated methods have been developed to 
solve this inverse problem in the framework of the single harmonic oscillator 
model [H El El El El El [101 [11] but they usually assume simple forms of the tip-surface 
force F ts in which the interaction depends on the instantaneous tip position only. 

2. The force disk 

2.1. General properties of surface forces 

Surface forces depending on the instantaneous tip position and velocity can be 
represented as two dimensional function in the z-z phase plane. The representation 
of a force depending on the history of the tip motion {z(t)}, for example hysteretic 
forces due to capillary formation or chemical bonds, is more difficult. However, for 
high quality factor oscillators with large stored energy, the motion is restricted to a 
narrow band near resonance and the steady state motion is well approximated by a 
sinusoidal trajectory. For a sinusoidal drive signal at the first resonance frequency 
the motion is then given by [23] 



where A is the oscillation amplitude, is the phase lag with respect to the drive. 
The motion orbits defined by equation (El) do not intersect each other in the z-z 




(4) 



z{t) = Acos(<4 >t + 4>) + h 



(5) 
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phase plane and thus every point in the phase plane can be mapped to a unique 
tip trajectory. We assume that the force along each trajectory does not depend on 
previous oscillation cycles. In this manner we can incorporate the dependence of the 
force on the past tip trajectory into the dependence on the instantaneous tip position 
and velocity. 

In every experiment there exists a maximum oscillation amplitude A max and 
velocity f max - We normalize the tip position and velocity to the maximum amplitude 
and velocity and subtract the static probe height, 

x = — (6) 

^max 

x = J- (7) 

to obtain a new force function on the closed unit disk in the x-x plane, 

F^ s \x,x) = F ts (A max x,v max x). (8) 

The motion defined by equation (jSJ) corresponds to circular orbits in the x-x plane 
with a maximum radius of 1 as shown in fig. [TJ Thus, it is sufficient to define the 

(x) 

force function F{ s as a function on the closed unit disk in the x-x plane. This force 
disk is a compact and general description of the interaction between the tip and the 
surface in narrow band dynamic AFM. 

One regularly distinguishes between conservative and non-conservative or 
dissipative surface forces. The force during a complete tip oscillation cycle can readily 
be decomposed into an effective conservative and an effective non-conservative force 
by finding the symmetric and anti-symmetric part of the measured force around the 
lower turning point of the tip motion[5l [7]. This decomposition ensures that the 
energy dissipation integral 

E dis = I F ts dz = f T F ts 0(f), z(t), {z(t)}) z(t) dt (9) 

JC[z(t)] JO 

equals zero for the effective conservative force. In contrast, dissipative forces are 
responsible for the energy dissipation. One should note that forces of different 
physical origin can contribute to the total surface force. Each of the effective forces 
may therefore contain contributions from forces which alone would be of purely 
conservative, purely dissipative or of mixed character. 
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Normalized position x 



Figure 1. Force disk in the x-x plane of a purely conservative force. For typical 
values of the static probe height and the oscillation amplitude the interaction is 
very localized at the surface. The depicted orbits correspond to purely sinusoidal 
tip motion. The orbits do not intersect each other so every point in the plane can 
be mapped to a unique oscillation amplitude. There is one value of Fj and Fq for 
each orbit which we assign to the lower turning point. 



In the same manner we are able to decompose the force disk into an effective 
conservative disk and into an effective dissipative disk Ffg), 

F t [ x \x,x) = F^(x,x) + F^(x,x), (10) 

by requiring that for the conservative force disk, no energy is dissipated during 
one complete oscillation cycle of the sinusoidal tip motion. This implies that the 
conservative force disk is symmetric with respect to the x-axis whereas the non- 
conservative force disk is anti-symmetric [10] 

F c (x,-x) =F c (x,x) (11) 
F BC (x,-x) = -F ac (x,x) (12) 
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Functions on the unit disk are naturally expressed in polar coordinates so we 
perform a change of variables to obtain a force function F t g that depends on the 
normalized amplitude r and the instantaneous phase of the oscillation 9, 

Fl x) (x,x)^Fg\r,0). (13) 

(r) 

Similar to functions on the unit sphere F ts can be expanded into a set of orthogonal 
functions such that 

oo n 

itW) = E E 4 m) 4 m) M)- (14) 

n=0 m=—n 

A natural choice of the basis functions are the Zernike polynomials |24j which 
can be defined as 

(r) cos(m9) for m > 

Zi m) = { n n " ( 15 ) 

[R^ ml) {r)sm{m9) for m < 

where the polynomials R^™ 1 ^ are given by[25j 

f (n "^ l)/2 (-mn-k)\ 

i2jM)( r ) = J ^ fc!(( n + m)/2-A;)!((n-m)/2-A;)! 

[ for 
and fulfill the orthogonality relation 

f R^\r)R^\r)r dr = ^—-J n ,n>- (17) 
Jo 2{n + 1) 

We will use this Zernike expansion of the force disk in the following section to 
investigate which parts of the force disk are measurable with narrow band dynamic 
AFM. 



2.2. Probing the force disk with dynamic narrow frequency band AFM 

In dynamic AFM the tip-surface force can also be considered as a time-dependent 
force acting on the oscillating tip. In AM-AFM and FM-AFM only the Fourier 
component of the time-dependent tip-surface force at the tip oscillation frequency 
is measurable above the noise floor. Higher frequency components of the force are 
filtered out by the high-quality factor resonance of the cantilever. The force Fourier 
component at the oscillation frequency can be expressed as a real-valued component 
Fj that is in-phase with a sinusoidal tip motion and a real-valued component Fq 
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that is quadrature to the tip motion[12]. With equation dSJ) we can write Fj and Fq 
as two integral equations 
1 

Fj = — J F ts (A cos (u§-k) + h, -utpAsin (w^t)) cos (u^k) d<18) 

Fq = T~Jo Fts ( A C ° S H 15 *) + ft ' H^)) sin H^) dt(19) 

The component Fi is the so-called virial of the tip motion which is only affected by 
the effective conservative force[23J. In contrast, Fq is connected to the dissipative 
interaction and comparison with the energy dissipation integral in equation (jHJ) yields 

E Ais = -2nAF Q . (20) 

Through their dependence on the tip motion z(t) and the tip velocity z(t) the force 
components Fj and Fq are functions of the oscillation amplitude A and the static 
probe height h. Alternatively, Ft and Fq can be considered as functions of h and 
the lower turning point z miQ such that at fixed static probe height, the force disk and 
the Fj(A) and Fq(A) curves share the same position axis as shown in figure [TJ 

While imaging with conventional dynamic AFM the feedback is working to 
keep the the oscillation amplitude constant. Thus, only one value for Fj and Fq 
is measured. Therefore, the combination of imaging and force measurement is not 
possible in conventional dynamic AFM and most force spectroscopy techniques rely 
on a measurement of the h dependence of Fj and Fq. As an alternative approach, we 
recently introduced the rapid measurement of the oscillation amplitude dependence of 
Fi and Fq at fixed static probe height with Intermodulation AFM[12J. To understand 
what information about the force disk can be extracted from a measurement of Fj(A) 
and Fq(A) we insert the Zernike expansion of the tip-surface force disk, equation ( fl4l) . 
into the integral equations (ITB1) and fTT9l 



Fj(A) = / J2 E at ) Rlr ) ( A /^)cos(me) C os(e)de 

Jo n=0 m=0 

oo 

^E^i^AVa*), (21) 

n=0 

2tt oo n 

Fq(A)= / E Y, ai n ) R { n\H^)sin{me)sin{e)de 

J Q 11=0 171=1 



tE^^H^). (22) 

n=0 
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We note that Fj(A) and Fq(A) only depend on the coefficients a^ 1 ^ with m = ±1 
where the individual coefficients can be recovered by using the orthogonality relation 
ffTTj) . This implies that the measurable information in narrow band dynamic AFM is 
fundamentally limited. It is not possible to reconstruct an arbitrary force disk from 
the experimentally available quantities Fj(A) and Fq(A). Nevertheless, the Fi(A) 
and Fq(A) curves provide qualitative insight into the interaction between tip and 
surface. 

The fundamental limitation at fixed static probe height is due to the limited 
frequency band. To increase the measurable information additional frequency bands 
have to be considered. One way to achieve this is to externally excite multiple 
cantilever eigenmodes |261 [271 128]. In liquid environments higher eigenmodes exhibit 
measurable response even in the absence of an external drive, due to the mode's low 
quality factor [291 130] . However, multiple eigenmodes require a higher dimensional 
phase space description of the dynamics and the representation of force must be 
appropriately adapted. Another way to introduce additional frequency components 
is the measurement of higher harmonics of the tip motion, although higher harmonics 
are only measurable under special conditions[6, 3 lj or with specialized cantilevers[32j. 



2.3. Examples of force disks 

In figure [2] four examples of the force disk for physically relevant model forces are 
shown in the z-z plane. Since the interaction is very localized we focus on to the 
region close to the surface. In figure a conservative van der Waals-Derjagin- 
Muller-Toropov (vdW-DMT) force is shown which is given by equation flA.ll) j 33] . 
The symmetry of the force disk with respect to the velocity and the fact that Fq 
equals zero for all motion turning points, is due to the vdW-DMT force being purely 
conservative. From the Fj(A) curve we can conclude that the conservative force can 
be decomposed in two regions of different character. Further away from the surface 
Fj is positive which corresponds to an average attractive force acting on the tip. 
When the tip makes contact with the surface Fj quickly becomes negative due to 
the rapid turn-on of the strongly repulsive force. 

Figure Eb shows a force disk that models a long-range dissipative force [34]. The 
model builds on a conservative vdw-DMT force with a Hamaker constant whose value 
is different for the tip approach and retract. This discontinuous behaviour is defined 
by equation ( 1A.2I) and appears in the force disk as a discontinuous asymmetry with 
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respect to the velocity. The force shown in figure [2Jd illustrates how the addition of 
a so-called dissipative interaction can modify the effective conservative tip-surface 
force. Despite the fact that the vdw-DMT parameters are the same as for figure 
Eh, the Fi curves differ for the two models. The added "dissipative" interaction 
actually contains an effective conservative component, which makes the total effective 
conservative force more attractive and less repulsive. The nonzero values of Fq 
indicate the presence of a dissipative force which is already present relatively far 
away from the surface. Moreover, from the shape of the Fq curve we can draw 
conclusions about the nature of the dissipative interaction: Further away from the 
surface Fq is inversely proportional to the squared oscillation amplitude. Below the 
contact point at z — 0, the derivative of Fq with respect to the oscillation amplitude 
reveals that Fq is inversely proportional to the oscillation amplitude in this region. 
This behaviour indicates that the energy dissipated during each oscillation cycle 
grows linearly with oscillation amplitude from the contact point on. 

In contrast to figure [2b the force displayed in figure [2b is an example of an 
additional dissipative interaction that does not modify the effective conservative 
force. The force is defined in equation flA.31) and combines a vdw-DMT force with a 
position-dependent viscous damping [35]. Due to the linear dependence of the force 
on the tip velocity the Fj curve in fig [2b does not differ from the Fj curve for the 
purely conservative vdW-DMT force in figure Eb . The damping coefficient for the 
force in figure [2b depends exponentially on position, which results in a dissipated 
energy that is proportional to the product of the oscillation amplitude and modified 
Bessel function of first kind in the oscillation amplitude [35J. Hence, the Fq curve 
also follows a modified Bessel function of first kind. 

Physically different interactions are encountered for capillary surface forces 
which can be due to adsorbed water layers on the surface. The common model is 
based on a vdW-DMT force with an additional hysteretic adhesion force that turns 
on and off instantaneously when the tip passes certain threshold positions as defined 
in equation (1A.4l) [36j. Up to a threshold position z on the force disk is symmetric 
with respect to velocity. At the threshold position the tip makes contact with the 
water layer on the surface which results in an additional attractive force. Due to 
this adhesion a capillary neck builds up between the tip and the surface when the 
tip retracts again. This necks breaks at a position z Q g that is different from z on . The 
water neck effectively extends the region of strong attraction during the tip retract. 
The abrupt force jumps are also visible in the Fi(A) and Fq(A) curves. Since the 
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dissipated energy is independent of the oscillation amplitude, the Fq(A) curve is 
inversely proportional to the oscillation amplitude which becomes more pronounced 
in the derivative of Fq with respect to the oscillation amplitude. 

For all example surface forces considered here the Fj(A) and Fq(A) curves 
give insight into the character of effective conservative and dissipative interactions 
between the tip and the surface. With ImAFM[37J [38] these curves are measured at 
each image point while scanning at imaging speeds[12]. This rapid and information- 
rich data acquisition technique allows for an enhanced interpretation of imaging 
contrast, a polynomial reconstruction of the tip-surface force [391 14*0] and the 
extraction of material properties |41j. 

3. Amplitude-dependence force spectroscopy (ADFS) 

The ultimate goal of dynamic AFM is the combination of high-speed and high- 
resolution imaging with high accuracy force measurements. However, imaging is 
performed at static probe height h above the surface and we have seen that in 
the case of fixed h it is not possible to reconstruct a complete force disk from 
the available Fj(A) and Fq(A) curves (see equations ( 1211) and ( l22l) ). Often we 
can make physically well-motivated assumptions about the tip-surface interaction 
which effectively correlate the Zernike expansion coefficients of the force disk. Under 
these assumptions a measurement of the amplitude dependence of Fj and Fq is 
sufficient for a quantitative reconstruction the tip-surface interaction. We call this 
method amplitude dependence force spectroscopy (ADFS) and in this section we 
show the reconstruction of the effective conservative tip-surface interaction as well 
as the reconstruction of a position-dependent viscous damping. 

3.1. Reconstruction of conservative tip-surface interactions 

In most cases it can be assumed that the effective conservative tip-surface force does 
not depend on tip velocity and tip motion history but is rather a function of the tip 
position only, 



We recently demonstrated how this assumption allows for the reconstruction of the 
effective conservative surface force F c with ADFS[42J. Using the equations (1241) and 




(23) 
(24) 
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Fj(A)= — F c (Acos9 + h)cos(9)d9 (25) 
2n Jo 

+ (26) 

Usually, the interaction length of the force is small (several nm) compared to the 

oscillation range (tens of nm). So we can neglect the upper half of the integration 
interval, 

F,(A)» f_ A F e(z + h)7 ^£= dz (27) 

•*<-^ + *> <h (28) 



2ttA 7^ -m 
With the definitions 

I = A 2 (29) 

F c (u) = F c (-^+h) (30) 

Fj(A)= -2n\fAF I (\fA) (31) 
equation (128]) becomes 

F 2 (l) = [ A 4== du - ( 32 ) 



We note that Fj(A) in equation ( 132]) is the Abel transform of the force F c (u). The 
Abel transform has a unique inverse [131 HI] with which we solve equation (1321 for 
the force F c (— z + /t) 

F c (u) = -— / —!LJ=dA 33 



i d r 2 VafaVa) ~ 

F c (-z + h)=-~— ^ J dA. 34 



In practice, the integral in equation ([34]) has to be evaluated numerically which is 
made difficult by the square root singularity at the upper integration limit. It is 
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therefore advantageous to perform the substitution y 2 = z 2 — A in the integral to 
remove the singularity and to improve the numerical stability, 

F c (-z) V ^ Z = A - -4- f Jz 2 - y 2 Fj (jz 2 - y 2 ] dy. (35) 



o 



z dz 

The resulting integral is then suitable for standard numerical integration methods 
like simple trapezoidal integration. 

3.2. Reconstruction of position- dependent damping 

Often dissipation is modeled as a position-dependent damping or friction coefficient 
X(z) such that the dissipative force is given by 

F ac {z,z) = \{z)z. (36) 

In this case we can reconstruct the effective damping function A (z) from the 
amplitude- dependence of Fq(A). We start with rewriting the integral equation ( TT9l) 
by using the equations (|T2|) and ( 136]) . 

uoA r 2 ^ 
2tt Jo 

z=Acos 9 w ' ' 
71 A J -A 

As for the conservative force we assume that the oscillation range is bigger than the 
interaction range of the dissipative force such that X(z) = for z > 0, 

F Q (A) « - ^- f° X(z + h)VA 2 - z 2 dz (39) 

TlA J -A 

"3" - f X{ -^ +h) ^—u du (40) 
Itx Ah s/u v ; 



Fn(A)= - ^— r \(Acos6 + h)sm 2 9 d9 (37) 

Z7T JO 

( A X(z + h)VA 2 -z 2 dz (38) 

J -A 



We use definition ff29|) and define additionally 

F Q( A^ - 2 ^A Fq (4A) (4!) 
A(„) s A( -^i +,i) (42) 



to arrive at 



A „ 







F Q (A)= X(u)^A-udu. (43) 
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The integral (j4"3j) represents a convolution of X(u) with y/u which can readily be 



solved in Laplace space (see Appendix B). For tip motion as defined by equation (jSJ) 



the position-dependent damping coefficient is then given by 
~ , ~, 2 d 2 f A F Q (u) , 

X(A) = --j- / -^LJ= du (44) 



X M A + h ) = -^^[ A ^m (45) 



3.3. Numerical results 

To validate the reconstruction of the effective conservative tip-surface force and the 
position dependent-damping we simulate the tip motion in a scalar force field with 
a vdW-DMT component and exponential damping as defined by equation ( 1A.3I) . 
To integrate the equation of motion (T4J) we use the CVODE solver with adaptive 
step-size and discrete event detection jl5] to account for the piecewise definition of 
the tip-surface force in equation (1A.3j) . We assume a cantilever typically used for 
experiments under ambient conditions with a quality factor of Q^ 1 ' = 400.0, a stiffness 
of fcW = 40.0 N/m and a resonance frequency of /q = 300.0 kHz. To probe the 
tip-surface interaction we use an ImAFM drive scheme such that the free oscillation 
amplitude is modulated between and 35 nm in a time window of T = 2 ms at a 
fixed static probe height of h = 22.0 nm above the sample surface as shown in figure 

Ek, 

From the tip motion close to the surface the Fj(A) and Fq(A) curves are 
obtained[12j which are then used for the numerical evaluation of the equations 
( |34|) and ( 1431) . The results are shown in figure |3j Both the reconstructed effective 
conservative force and the position-dependent damping are in excellent agreement 
with the force and damping actually used in the simulation. For the conservative force 
the actual and the reconstructed force deviate by at most 4.5 % of the maximum force. 
This deviation is the result of the rapid change in the force around the contact point, 
which generate Fourier components with amplitudes typically below the noise level 
in AFM systems and these components have been neglected in the reconstruction. 

We note that the interaction is localized in the region between -3.2 nm and 
7.5 nm. A more optimized drive scheme would therefore modulate the oscillation 
amplitude only in this interaction region, thereby sampling the interaction more 
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smoothly since more turning points lie in the interaction region during the same 
measurement time [1H1 H7] . 

4. Conclusions 

We have introduced a comprehensive framework for the representation of tip-surface 
interaction at fixed static probe height in terms of a "force disk". This framework 
allows for a complete description of forces depending on the instantaneous tip 
position and velocity, appropriate to narrow band dynamic AFM. We show that 
it not possible to fully reconstruct arbitrary tip-surface interactions from narrow 
band dynamic AFM measurements. However, by considering different tip oscillation 
amplitudes, the experimental curves Fj(A) and Fq(A) give insight in the character 
of the effective conservative and the effective dissipative interaction. With minimal 
additional assumptions these curves can be used for quantitative reconstructions 
of the effective conservative tip-surface force and the effective position-dependent 
damping. 

The framework introduced here gives a more solid theoretical foundation for 
the visualization and analysis of tip-surface forces than traditional force curves. 
Furthermore, we expect that the notion of a force disk will inspire new multifrequency 
force probing schemes expanding the capabilities of conventional narrow band 
dynamic AFM. On the experimental side amplitude-dependence force spectroscopy 
enables rapid force measurements while scanning sample surfaces at normal imaging 
speeds. From the obtained force volume data sets different mechanical properties of 
the sample surface can be derived without the constraints and limitations of over- 
simplified interaction models. 

Since the underlying physics of the theory presented here is the classical 
harmonic oscillator moving in an external force field. We see wide-spread applications 
of this method of analysis in other fields using resonant detection techniques based 
on different kinds of high-quality factor resonators [4"6l Wf\ . 
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Appendix A. Definition of model surface forces 



In section 12.31 we introduced various force models of which the mathematical 
definitions are given here. Also given are the values of the parameters used when 
plotting or simulating with these models. 

The basis of all models forms a vdW-DMT force which is defined as 

HR 



Fdmt(z) 



6(a 
HR 



for z > 



4 



+ -E*y/R(-z) 3/2 for z < 



(A.l) 



Gal ' 3 

where H = 3.28 • 1CT 17 J is the Hamker constant, R = 10.0 nm is the tip radius, 
ao = 2.7 nm is the intermolecular distance and E* = 1.5 GPa is the effective stiffness 
of the tip-surface system. 

The long-range dissipative force is defined as 



HappR 



6(a - z) 2 
6ao 3 

H ret R 



6(a - z) 



6<2n ' 3" 



for z > and z < 
for z < and z < 
for z > and z > 
for z < and z > 



(A.2) 



which extends the vdW-DMT model by a Hamaker constant that depends on if the 

tip is approaching the surface (H app = H) or is retracting from it (H Tet = 2H) 

The vdW-DMT with an additional viscous damping term depending 

exponentially on position is given by 

HR . . . 

7o exp {— z /z 1 ) z for z > 



Ee^pi^Z, Z^j 



6(a - zf 
HR 4 

— -J- + -^v^-z) 372 - 7o exp (-*/*,) i for z < 
6ao 3 



(A.3) 



where 70 = 3.5 Js is the damping factor and z 7 
length. 



1.5 nm is the damping decay 
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Capillary interactions are modeled as 

HR 



F cap (z,{z(t)}) 



6(a 



HR 



6(a 



HR 



6(a - zf 
HR 



47T7h 2 o^ 
1 + 



^H 2 

47T7h 2 o-R 



6(a - zf 1 + 



/lH 2 



n3/2 



1 + 



_ao 

/»H 2 



for z > z s 

for z > z on and z < z s an( 
for z > z on and z < z g an< 

for z < z on and z > 

for z < 



where 7h 2 o = 72 • 10~ 3 J/m 2 is the surface energy of water and hn 2 o = 0.1 nm is the 
effective thickness of the adsorbed water layer. The state variable m is set to 1 when 
the tip makes contact with the water layer at z on = 0.2 nm and is set to when 
the water neck breaks at z q q = 1.94 nm. In contrast to the model introduced injSS] 
we assume that the interface between the surface and the water layer is located at 
z = nm. 

Appendix B. Inversion of the damping integral 

The reconstruction of the position-dependent damping X(z) from Fq(A) requires in 
the inversion of the integral equation 

F Q (A) = f \{u)^A-udu (B.l) 
Jo 

We study this convolution in Laplace space where it becomes 

C {f q (A)} = C ^J* \{u)^A-u du\ = C {A(I)} C { ^3} . (B.2) 

The Laplace transform of the square root is £ = ^s~ 3 ^ 2 which yields 

^C{\(A)} =h-^c{F Q (A)} (B.3) 

= I c {ti} c ^ q(A) } (R4) 
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(B.5) 



Now, we use the relation between integration in the time domain and its Laplace 
transform and obtain 
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Figure 2. Model force disks in the z-z plane at a static probe height of h = 22 nm 
and a maximum oscillation amplitude A max = 25 nm. The models all build 
on a purely conservative vdW-DMT force (a) and add long range adhesion(b), 
position-dependent viscous damping(c) and capillary interactions(d) respectively. 
We concentrate on the region close to the surface and display the disk sections 
together with the corresponding Fj and Fq curves as functions of the lower motion 
turning point. 



Tip-surface interactions in dynamic atomic force microscopy 



21 




5 10 15 
Position z [nm 



5 10 15 20 
Position z [nm] 



Figure 3. Simulated tip motion in the time domain(a). On a fast time scale the tip 
motion is nearly purely sinusoidal with an amplitude that is modulated on a much 
slower time scale. The ADFS reconstructions of the effective conservative force 
curve(a) and the position-dependent damping curve(b) are in excellent agreement 
with the reference curves used in the simulation. 



